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We investigate heating of the cool core of a galaxy cluster through the dissi- 
pation of sound waves and weak shocks excited by the activities of the central 
active galactic nucleus ( AGN) . Using a weak shock theory, we show that this 
heating mechanism alone cannot reproduce observed temperature and den- 
sity profiles of a cluster, because the dissipation length of the waves is much 
smaller than the size of the core and thus the wave energy is not distributed 
to the whole core. 

1 Introduction 

The failure of standard cooling flow models indicates that the gas is prevented 
from cooling by some heating sources. At present, the most popular candi- 
date for the heating source is the active galactic nucleus (AGN) at the cluster 
center. However, it is not understood how the energy ejected by the AGN 
is transfered into the surrounding ICM. One idea is that bubbles inflated by 
AGN jets move outward in a cluster by buoyancy and mix the surrounding 
ICM [1, 2, 3]. As a result of the mixing, hot ICM in the outer region of the 
cluster is brought into, and subsequently heats, the cluster center. The other 
idea is that the dissipation of sound waves created through the AGN activi- 
ties. In fact, sound waves or weak shocks that may have evolved from sound 
waves are observed in the Perseus and the Virgo clusters [4, 5] . It was argued 
that the viscous dissipation of the sound waves is responsible for the heating 
of a cool core [4, 6]. They estimated the dissipation rate assuming that the 
waves are linear. However, when the amplitude of sound waves is large, the 
waves rapidly evolve into non- linear weak shocks [7] , and their dissipation can 
be faster than the viscous dissipation of linear waves. Although it was argued 
the presence of weak shocks in Ref. [4], their evolution from sound waves 
was not considered. Numerical simulations of dissipation of sound waves cre- 
ated by AGN activities were also performed [8] . Their results actually showed 
that the sound waves became weak shocks. However, their simulations were 
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finished before radiative cooling became effective. Thus, the long-term bal- 
ance between heating and cooling is still unknown. In this paper, we consider 
the evolution of sound waves to weak shocks, and analytically estimate the 
'time-averaged' energy flux of the propagating waves as a function of dis- 
tance, explicitly taking account of the dissipation at weak shock fronts and 
its global balance with radiative cooling. We assume the Hubble constant of 
Hq = 70 km s _1 Mpc -1 . The detail of the models and results are shown in 
Ref. [9]. 



2 Models 

We assume that sound waves are created by central AGN activities. The 
waves propagate in the ICM outwards. These waves, having a relatively large 
but finite amplitude, eventually form shocks to shape sawtooth waves. If the 
velocity amplitude is larger than ~ 0.1 sound velocity (the Mach number is 
^ 1.1), those waves steepen and become weak shocks after propagating less 
than a few wavelengths [10]. These shock waves directly heat the surrounding 
ICM by dissipating their wave energy. We adopt a heating model for the solar 
corona based on a weak shock theory ([10, 7], see also [11]). We assume that 
a cluster is spherically symmetric and steady. 
The equation of continuity is 

M = -4Trr 2 pv , (1) 

where M is the mass accretion rate, r is the distance from the cluster center, 
p is the ICM density, and v is the ICM velocity. The equation of momentum 
conservation is 

dv = GM(r) I dp 1 1 d 2 

dr ' r 2 pdr pc s {l + [(7 + l)/2]a w } r 2 dr 1 w> K> 

where G is the gravitational constant, M{r) is the mass within radius r, 
p is the ICM pressure, c s is the sound velocity, j(= 5/3) is the adiabatic 
constant, and a w is the wave velocity amplitude normalized by the ambient 
sound velocity (a w = Sv w /c s ). The wave energy flux, F w , is given by 

1 .3 2 ( - T+l 



Fw = ^pcia z w ^1 + ~^-a w j . (3) 

The energy equation is 

d (\ 2 7 k B T\ GM(r) 1 d r 2/rn „ N1 , „ 

pv-r ^v 2 + -^-^^ + P v 5^ + - — [r 2 (F w + F c )]+n 2 e A(T) = , 

dr \ 2 7 — 1 prriH ) r z r z dr 

(4) 

where ks is the Boltzmann constant, T is the ICM temperature, p(= 0.6) is 
the mean molecular weight, ran is the hydrogen mass, F c is the conductive 
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flux, n e is the electron number density, and A is the cooling function. The term 
V • F w indicates the heating by the dissipation of the waves. The equation 
for the time-averaged amplitude of the shock waves is given by 



da w 
dr 



Idp _ 2(7+1)0^ _ 2 
p dr c s t r 



1 dc s 
c, dr 



(5) 



where r is the period of waves, which we assume to be constant [7, 10]. The 
second term of the right side of equation (5) denotes dissipation at each shock 
front. 



3 Results 

For parameters of our model cluster, we adopt the observational data of the 
Perseus cluster [12] . We assume that r s = 280 kpc, M(r woo ) = 3.39x 1O 14 M , 
and riooo = 826 kpc, where the mean density within r^ is A times the 
critical density of the Universe. Waves are injected at the inner boundary 
r = ro, which should be close to the size of bubbles observed at cluster 
centers. We assume that Ao = ro, where Ao is the initial wavelength. If the 
waves are injected in a form of sound waves with amplitude 0.1 ^ a w < 1, 
waves travel about Ao before they become shock waves [10]. Therefore, for 
r < r < r + A = 2Ao, we assume that V • F w — (eqs. [2] and [4]), 
and that the second term of the right-hand side of equation (5) is zero. The 
temperature, electron density, and wave amplitude at r = r are T n e0 , 
and a W Q, respectively. Unless otherwise mentioned, the first two are fixed at 
T = 3 kcV and n e0 = 0.08 cm~ 3 , respectively, based on the observational 
results of the Perseus cluster [13]. 
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Fig. 1. (a) Temperature, (b) density, and (c) dissipation length profiles for a w o = 
0.3 (dotted lines), 0.5 (thin solid lines), and 0.7 (dashed lines). Other parameters are 
r = 1 x 10 7 yr, M = 5OM0yr _1 , and f c = 0. Filled circles are Chandra observations 
of the Perseus cluster [13]. The bold solid lines correspond to a genuine cooling flow 
model of M = 500 M Q yr~\ 
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In Figure 1, we show the results when r is fixed at 1 x 10 7 yr, and a w0 
is changed. For the Perseus cluster, it was estimated that a w o ~ 0.5 [4]. The 
dissipation length is defined as l w — \F W /V ■ F w \. For these parameters, the 
initial wavelength is Ao = 9 kpc, which is roughly consistent with the Chandra 
observations [4]. Other parameters are M — 50 M Q yr^ 1 , and f c — 0, where 
f c is the ratio of actual thermal conductivity to the classical Spitzer con- 
ductivity. In general, larger M reproduces observed temperature and density 
profiles better. However, large M is inconsistent with recent X-ray obser- 
vations as was mentioned in § 1. For comparison, we show the results of a 
genuine cooling flow model (M = 500M Q yr~ 1 , a w0 — 0, and f c — 0) and the 
Chandra observations of the Perseus cluster [13]. Figures la and lb show that 
only a small region is heated. The jumps of T and n e at r = 2Ao = 18 kpc are 
produced by weak shock waves that start to dissipate there. The energy of the 
sound waves rapidly dissipates at the shocks, which is clearly illustrated in 
short dissipation lengths, l w ~ 2-15 kpc (Fig. lc). These dissipation lengths 
are smaller than those of viscous dissipation for linear waves, which can be 
represented by l v = 420 XqUq.q^T^ 2 kpc, where the wavelength A = 9Ag kpc, 
the density n = 0.08 no.os cm~ 3 , and the temperature T = 3 T 3 keV [4]. In 
Figure 1, the ICM density becomes large and the temperature becomes small 
at r ^ 2Ao so that the rapid shock dissipation is balanced with radiative 
cooling. Because of this, waves cannot reproduce the observed temperature 
and density profiles that gradually change on a scale of ~ 100 kpc. 
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Fig. 2. Same as Fig. 1 but for r = 1 x 10 7 yr (solid lines), and 2 x 10 7 yr (dashed 
lines). Other parameters are a w o = 0.5, M = 50 M yr -1 , and f c = 0. 



In Figure 2, we present the results when t = 2 x 10 7 yr. Compared with the 
case of t = 1 x 10 7 yr, the wave energy dissipates in outer regions. However, 
the dissipation lengths are still smaller than the cluster core size (~ 100 kpc). 
Note that larger r (or Ao) means formation of larger bubbles. As indicated 
in Ref. [14], it is unlikely that the size of the bubbles becomes much larger 
than 20 kpc; the bubbles start rising through buoyancy before they become 
larger. On the other hand, when r < 10 7 yr (A < 9 kpc), the waves heat only 
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the ICM around the cluster center. The predicted temperature and density 
profiles are obviously inconsistent with the observations. 

The inclusion of thermal conduction changes the situation dramatically. 
Figure 3 shows the results when f c = 0.2. The models including both wave 
heating and thermal conduction can well reproduce the observed temperature 
and density profiles. Figure 3c shows the contribution of the wave heating 
(— V • F w ) to compensating radiative cooling {n 2 e A). Since —V • F w /n 2 e A > 
1/2 for r ~ 20-30 kpc, the wave heating is more effective than the thermal 
conduction in that region. 
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Fig. 3. (a) Temperature, (b) density, and (c) dissipation strength profiles for a w o = 
and To = 3 keV (dotted lines), a w o = 0.5 and To = 3 keV (solid lines), and 
a w o = 0.5 and To = 2.5 keV (dashed lines). Other parameters are r = 1 x 10 7 yr, 
M — 50 Mq yr -1 , and f c = 0.2. Filled circles are Chandra observations of the 
Perseus cluster [13]. 



4 Discussion 

We showed that sound waves created by the central AGN alone cannot re- 
produce the observed temperature and density profiles of a cluster, because 
the dissipation length of the waves is much smaller than the size of a clus- 
ter core and the waves cannot heat the whole core. The results have been 
confirmed by numerical simulations [15]. On the other hand, we found that 
if we include thermal conduction from the hot outer layer of a cluster with 
the conductivity of 20% of the Spitzer value, the observed temperature and 
density profiles can be reproduced. The idea of the "double heating" (AGN 
plus thermal conduction) was proposed in Ref. [16]. 

However, the fine structures observed in cluster cores may show that the 
actual conductivity is much smaller than that we assumed [17]; the structures 
would soon be erased, if the conductivity is that large. If the conductivity is 
small, we need to consider other possibilities. While we considered successive 
minor AGN activities, some authors consider that rare major AGN activities 
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should be responsible for heating of cool cores [18]. In this scenario, powerful 
bursts of the central AGN excite strong shocks and heat the surrounding gas 
in the inner region of a cluster on a timescale of ^ 10 9 yr. Moreover, in this 
scenario, heating and cooling arc not necessarily balanced at a given time, 
although they must be balanced on a very long-term average. This is consis- 
tent with the fact that there is no correlation between the masses of black 
holes in the central AGNs and the X-ray luminosities of the central regions 
of the clusters [19]. Alternative idea is that cluster mergers are responsible 
for heating of cool cores [20, 21]. In this "tsunami" model, bulk gas motions 
excited by cluster mergers produce turbulence in and around a core, because 
the cooling heavy core cannot be moved by the bulk gas motions, and the 
resultant relative gas motion between the core and the surrounding gas in- 
duces hydrodynamic instabilities. The core is heated by turbulent diffusion 
from the hot outer region of the cluster. Since the turbulence is produced and 
the heating is effective only when the core is cooling and dense, fine-tuning 
of balance between cooling and heating is alleviated for this model. 
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